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Abstract: 

We consider a linearly elastic composite medium, which consists of a homogeneous matrix containing statistically 
inhomogeneous random set of heterogeneities and loaded by inhomogeneous remote loading. The new general 
integral equation is obtained by a centering procedure without any auxiliary assumptions such as, e.g., effective 
field hypothesis implicitly exploited in the known centering methods. The method makes it possible to abandon 
the basic concepts of classical micromechanics such as effective field hypothesis, and the hypothesis of "ellip- 
soidal symmetry" . The results of this abandonment leads to detection of some fundamentally new effects that is 
impossible in the framework of a classical background of micromechanics. 
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1 Introduction 

The final goals of micromechanical research of composites involved in a prediction of both 
the overall effective properties and statistical moments of stress-strain fields are based on the ap- 
proximate solution of exact initial general integral equations connecting the random stress fields 
at the point being considered and the surrounding points. These equations are well-known for 
statistically homogeneous composite materials subjected to homogeneous boundary conditions. 
In the current paper, these known equations are generalized to the case of inhomogeneity of both 
the statistical microstructure and applied loading. The method is based on a centering proce- 
dure of subtraction from both sides of a known initial integral equation the statistical averages 
obtained without any auxiliary assumptions such as, e.g., effective field hypothesis implicitly 
exploited in the known centering methods. 

A considerable number of methods are known in the linear elasticity theory of composites. 
Appropriate, but by no means exhaustive, references are provided by the book by Buryachenko 
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[1]. Nowadays, it appears that variants of the effective medium by Kroner [2] and mean field 
methods by Mori and Tanaka [3] are the most popular and widely used methods. Recently, a 
new method has become known, the multiparticle effective field method (MEFM) (see [1]). The 
MEFM is based on the theory of functions of random variables and Greens functions. Within this 
method one constructs a hierarchy of statistical moment equations for conditional averages of 
the stresses in the inclusions. The hierarchy is then cut by introducing the notion of an effective 
field. This way the interaction of different inclusions is taken into account. Buryachenko [1] has 
demonstrated that the MEFM includes as particular cases the well-known methods of mechanics 
of strongly heterogeneous media. 

However, all mentioned methods are based on the effective field hypothesis (EFH, even 
if the term 11 effective field hypothesis" was not indicated) according to which each inclusion 
is located inside a homogeneous socalled effective field (see for references [1]). Effective field 
hypothesis is apparently the most fundamental, most prospective, and most exploited concept 
of micromechanics. This concept has directed a development of micromechanics over the last 
sixty years and made a contribution to their progress incompatible with any another concept. 
The idea of effective field dating back to Mossotti [4] was added by the hypothesis of "ellipsoidal 
symmetry" 1 for the distribution of inclusions attributed to Willis [5]. However, we will show in 
this paper that the EFH (also called the hypothesis HI) is a central one and other concept plays 
a satellite role providing the conditions for application of the EFH. Moreover, we will show that 
all mentioned hypotheses are not really necessary and can be relaxed. 

The outline of the study is as follow. In Section 2 we present the basic equation of thermoe- 
lasticity, notations, and statistical description of the composite microstructure. The new general 
integral equations are proposed in Section 3 for the case of statistically inhomogeneous structures 
of composite materials. These equations are obtained by a centering procedure of subtraction 
from both sides of a known initial integral equation the statistical averages obtained without any 
auxiliary assumptions such as, e.g., EFH implicitly exploited in the known centering methods. 
The new general integral equation is compared with the known ones. In Section 4 we recall the 
basic concepts defining the background of classical micromechanics. Explicit formulae for both 
effective elastic moduli and strain concentrator factor are presented. The new general integral 
equations presented in Section 5 through the operator form of the particular solutions for one 
heterogeneous in the infinite matrix subjected to inhomogeneous effective field. This equation 
is solved by the iteration method in the framework of the quasi-crystallite approximation but 
without basic hypotheses of classical micromechanics such as both the EFH and "ellipsoidal 
symmetry" assumption. In Section 5 we qualitatively explain the advantages of the new ap- 
proach with respect to the classic ones and demonstrate the corrections of popular propositions 
obtained in the framework of the old background of micromechanics. 

2 Description of the mechanical properties and geometrical structure of the com- 
ponents 

Let a linear elastic body occupy an open bounded domain w C R d with a smooth boundary 
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r and with a characteristic function W and space dimensionality d [d = 2 and d = 3 for 2- 
D and 3-1? problems, respectively). The domain w contains a homogeneous matrix and a 
statistically inhomogeneous set X = Vi,Xj) of inclusions t>j with characteristic functions Vi 
and centers Xj. It is assumed that the inclusions can be grouped into component (phase) 
with identical mechanical and geometrical properties (such as the shape, size, orientation, and 
microstructure of inclusions). For the sake of definiteness, in the 2-D case we will consider a 
plane-strain problem. At first no restrictions are imposed on the elastic symmetry of the phases 
or on the geometry of the inclusions t . 

We will consider the local basic equations of elastostatics of composites 

Vo-(x) = 0, (2.1) 
<r(x) = L(x)e(x), or e(x) = M(x)<r(x), (2.2) 
e(x) = [V <g> u + (V <g> u) T ]/2, Ince(x) = 0, (2.3) 

where <8> denotes tensor product, and (.) T denotes matrix transposition. L(x) and M(x) = L(x) _1 
are the known stiffness and compliance fourth-order tensors, and the common notation for con- 
tracted products has been employed. 

In particular, for isotropic constituents the stiffness tensor L is given in terms of the local 
bulk modulus k and the shear modulus fi: L = (dk, 2fi) = dfcNi + 2/xN2, Ni = 8 <8> S/d, N2 = 
I — Ni (d = 2 or 3); S and I are the unit second-order and fourth-order tensors. The tensors g 
(g = L, M) of material properties are piecewise constant and decomposed as g = g(°) + gi(x) = 
g(0) _|_ g^(x) where g( ) =const, g(x) = g(°) at x € t)° and g^(x) = g^ is a homogeneous 
function of the x G v^: 

L^(x) = = const, at x G . (2.4) 

The upper index of the material properties tensor put in parentheses shows the number of the 
respective constituent. The upper index ( m ) indicates the components and the lower index i 
indicates the individual inclusions; = w\v, v = V(x) = = X)^( x )> an d V^)(x) 
and Vi(x) are the indicator functions of and Vi, respectively. 

The boundary conditions at the interface boundary will be considered together with the 
mixed boundary conditions on T with the unit outward normal n r 

u(x) = u r (x), x G T u , (2.5) 
<r(x)n r (x) = t r (x), x G T t , (2.6) 

where T u and Tt are prescribed displacement and traction boundaries such that T u U Tt = 
T, T u nT t = 0. u r (x) and t r (x) are, respectively, prescribed displacement on T u and traction 
on iy Of special practical interest are the homogeneous boundary conditions 

u r (x) = e r x, e r = const., x G T, (2.7) 
t r (x) = cr r n r (x), cr r = const., x G T, (2.8) 

^It is known that for 2-D problems the plane-strain state is only possible for material symmetry no lower than 
orthotropic (see e.g. [6]) that will be assumed hereafter in 2-D case. 
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where e r (x) = ^[V <g> u r (x) + (V <8> u r (x)) T ], x G T, and e r and <x r are the macroscopic strain 
and stress tensors, i.e. the given constant symmetric tensors. 

It is assumed that the representative macrodomain w contains a statistically large number 
of realizations a of inclusions »j E » (providing validity of the standard probability theory 
technique) of the constituent Vi G v (i = 1, 2, . . .). A random parameter a belongs to a sample 
space A, over which a probability density p(a) is defined (see, e.g., [7]). For any given a, any 
random function g(x, a) (e.g., g = V,cr,e) is defined explicitly as one particular member, with 
label a, of an ensample realization. Then, the mean, or ensemble average is defined by the angle 
brackets enclosing the quantity g 

(g)(x) = / g(x,a)p{a)da. (2.9) 
J A 

No confusion will arise below in notation of the random quantity g(x, a) if the label a is dropped 
for compactness of expressions unless such indication is necessary. One treats two material length 
scales (see, e.g., [8]): the macroscopic scale L, characterizing the extent of w, and the microscopic 
scale a, related with the heterogeneities V{. Moreover, one supposes that applied field varies on 
a characteristic length scale A. The limit of our interests for both the material scales and field 
one is presented in an asymptotic sense 

L > A > a, (2.10) 

as the scale of microstructure a relative to the macroscale L tends to zero. All the random 
quantities under discussion are described by statistically inhomogeneous random fields. For the 
alternative description of the random structure of a composite material let us introduce a con- 
ditional probability density (p(yi, Xj|ui, xi, . . . , v n , x n ), which is a probability density to find the 
i-th inclusion with the center Xj in the domain Vi with fixed inclusions v\ , . . . , v n with the 
centers xi, . . . ,x n . The notation <p(vi, Xj|; v±, xi, . . . , v n ,x n ) denotes the case Xj ^ xi, . . . ,x n . 
We will consider a general case of statistically inhomogeneous media with the homogeneous 
matrix (for example for so-called Functionally Graded Materials (FGM)), when the condi- 
tional probability density is not invariant with respect to translation: (p(vi,Xi\vi,x±, . . . ,v n ,x n ) 
7^ ip(vi,Xi + x\v±, xi + x, . . . , v n , x n + x), i.e. the microstructure functions depend upon their ab- 
solute positions. In particular, a random medium is called statistically homogeneous in a narrow 
sense if its multi-point statistical moments of any order are shift-invariant functions of spatial 
variables. Of course, <p(vi, Xj|; v\, xi, . . . , v n ,x n ) = for values of Xj lying inside the "excluded 
volumes" Uf^j (since inclusions cannot overlap, m = 1, . . . ,n), where v^i D v m with character- 
istic function is the "excluded volumes" of x« with respect to v m (it is usually assumed that 
v mi = v m), and <p(y i ,x i \;v 1 ,x 1 , . . .,v n ,x n ) -> </?(^,Xj) as |xj - x m I -> oo, m = 1, . . . ,n (since 
no long-range order is assumed). (p(vi,x) is a number density, n = n(x) of component v 3 v i 
at the point x and c^\x) is the concentration, i.e. volume fraction, of the component Vi £ v 
at the point x: cW(x) = {V}(x) =U i n( 1 )(x), v { = mesvf, i = 1,2,...), c(°)(x) = 1 - (V)(x). 
Hereafter only if the pair distribution function g(xi — x m ) = Lp{vi,Xi\;v m ,x m )/n^ kS} depends 
on |x m — Xj| it is called the radial distribution function (RDF). The notations ((.))(x) and 
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((.)|t>i,xi; . . . ;v n ,x n )(x) will be used for the average and for the conditional average taken for 
the ensemble of a statistically inhomogeneous field X = {vi) at the point x, on the condition 
that there are inclusions at the points xi, . . . ,x n and Xj 7^ Xj if i ^ j = 1, . . . , n). The 
notations ((.)|; ui,xi; . . . ;w n ,x n )(x) are used for the case x ^ v±, . . . , v n . The notation ((-))j( x ) 
at x G Vi means the average over an ensemble realization of surrounding inclusions (but not over 
the volume V{ of a particular inhomogeneity, in contrast to ((•))(«)) at the fixed V{. Without loss of 

generality, we assume that the subdomains ifi^ do not touch the boundary T; such subdomains 
are called floating subdomains. In other words, the body w is considered as one cut out from an 
infinite random medium and the inclusions Vi intersected with the boundary V are replaced by 
the matrix material. 

We will use two sorts of conditional averages of some tensor g (e.g., g = e, cr). At first, the 
conditional statistical average in the inclusion phase (g)^(x) = (gV)^ q \yi) (at the condition 
that the point x is located in the inclusion phase x G v^) can be found as (gV)^ q \x.) = 
(^^(x)} -1 (gV^) (x) . Usually, it is simpler to estimate the second conditional averages of these 
tensors in the concrete point x of the fixed inclusion x € v q : (g|v 9 ,x 9 )(x) = (g) (x). At first 
we built some auxiliary set v q (x) with the boundary dv q (x) formed by the centers of translated 
ellipsoids v q (0) around the fixed point x. We construct v q (x) as a limit v® q — > i^(x) if a fixed 
ellipsoid Vk is shrinking to the point x. Then we can get a relation between the mentioned 
averages [x = (x 1 , . . . , x d ) T }: 

(g)^(x)= / ^)(y)(g|^,y}(x)dy. (2.11) 

Jv\ (x) 

Formula (2.11) is valid for any material inhomogeneity of inclusions of any concentration in the 
macrodomain w of any shape (if w^(x) C w). Obviously, the general Eq. (2.11) is reduced to 

the popular one (g)^ = (g) q for statistically homogeneous media subjected to homogeneous 
boundary conditions. 

3 General integral equation 

Substituting the constitutive equation (2.1) and the Cauchy equation (2.3) into the equilib- 
rium equation (2.1), we obtain a differential equation with respect to the displacement u which 
can be reduced to a symmetrized integral form after integrating by parts (see, e.g, Chapter 7 in 
[1]) 

e(x) = e°(x)+ / U(x-y)T(x)dy, (3.1) 

J w 

where r(x) = Li(y)e(y) is called the stress polarization tensor, and the surface integral is absent 
because the heterogeneous are assumed (without loss of generality) to be floating ones. The 
integral operator kernel U is an even homogeneous a generalized function of degree — d defined 
by the second derivative of the Green tensor G: C/y^(x) = [VjV/Gj^x)]^^, the parentheses in 
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indices mean symmetrization. G is the infinite-homogeneous-body Green's function of the Navier 
equation with an elastic modulus L (0) defined by V [l (0) V <g> G (0) (x) + (V G (0) (x)) T J /2 j = 
—65 (it), and vanishing at infinity (|x| — > oo), <5(x) is the Dirac delta function. The deterministic 
function e °(x) is the strain field which would exist in the medium with homogeneous properties 
and appropriate boundary conditions (see, e.g, [9]): 

4(x) = jf [G i(P) ,)(x-s)*?(s) - ^(s)4°L4°(L)/( x " s)n,-(s)]ds, (3.2) 

which conforms with the stress field <x°(x) = l/°)e°(x). The representation (3.2) is valid for 
both the general cases of the first and second boundary value problems as well as for the 
mixed boundary- value problem (see for references [1]). For simplicity we will consider only 
internal points x G w of the microinhomogeneous macrodomain w at sufficient distance from 
the boundary 

a <C |x-s|, Vs € T. (3.3) 

In so doing, some Cauchy data [u°(s), t°(s)] (3.2) (if they are not prescribed by the bound- 
ary conditions) will depend on perturbations introduced by all inhomogeneities, and, therefore 
e°(x) = e°(x,a). 

Now we will center Eq. (3.1), i.e. from both sides of Eq. (3.1) their statistical averages are 
subtracted 

e(x) = (s)(x) + f [U(x - y)r(y)-(U(x - y)r>(y)]dy + X\. (3.4) 

J w 

Without loss of generality, it is assumed the traction boundary conditions (2.6). Then = 
e (x, a) — (e )(x) defined by the integral (3.2) (containing only u (s)) over the external surface 
s G r can be dropped out, because this tensor vanishes at sufficient distance x from the boundary 
T (3.3) (see for details, e.g., Ref. [10] and its applications Shermergor [11]; see also [1], [12] where 
the case of nonfloating subdomains is considered). 

The integrals in Eqs. (3.4) converges absolutely for both the statistically homogeneous and 
inhomogeneous random fields X of inhomogeneities. Indeed, even for the FGMs, the term in the 
square brackets in Eq. (3.4) is of order 0(|x — y\~ 2d ) as |x — y| — > oo, and the integral in Eq. 
(3.4) converges absolutely. Therefore, for x G w considered in Eqs. (3.4) and removed far enough 
from the boundary T (3.3), the right-hand side integrals in Eq. (3.4) does not depend on the 
shape and size of the domain w, and it can be replaced by the integrals over the whole space 
R d . With this assumption we hereafter omit explicitly denoting R d as the integration domain in 
the equation 

e(x) = (s)(x)+ J [U(x-y)T(y)-(U(x-y)r(y))(y)]dy. (3.5) 
It should be mentioned that a popular equality 



(U(x-y)r(y))(y) = U(x-y)(r(y))(y) 



(3.6) 
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is only asymptotically valid at |x — y| — > oo. Then Eq. (3.5) is asymptotically reduced to the 
known one (see for details [1]) 

e(x) = (s)(x) + |u(x - y)[r(y) - (r)(y)]dy, (3.7) 

which in turn coincides with the equation 

e(x) = <£)+/ U(x - y)[ T (y) - (r)]dy, (3.8) 

for statistically homogeneous media subjected to the homogeneous boundary conditions. 

Let the inclusions v±,. . . ,v n be fixed and we define two sorts of effective fields £j(x) and 
£i,...,n( x ) {i = 1, • • • ,n; x G t>i, . . . ,u n ) by the use of the rearrangement of Eq. (3.6) in the 
following form (see for the earliest references of related manipulations [1]): 



e(x) = ei(x) + J U(x - y)Vi(y)T(y)dy, 
5i(x) = ei,...,„(x) + X; / U(x-y)^(y)r(y)dy, 



£i,...,n(x) = (e)(x) + J |u(x-y)r(y)F(y|;ui,xi;...;v n ,x n ) 

-(U(x-y)r)(y)}dy, (3.9) 

for x G fj, i = 1,2, ...,n; here V(y|ui, Xi; . . . ; v n , x n ) is a random characteristic function of 
inclusions x € w under the condition that Xj 7^ Xj ii i ^ j = 1, . . . , n). Then, considering 
some conditional statistical averages of the general integral equation (3.5) leads to an infinite 
system of new integral equations (n = 1, 2, . . .) 

n „ 

(e|«i,xi; . . . ;u n ,x n )(x) - ^ / U(x - y)(Fi(y)r|ui, xi; . . . ; u n , x„)(y)dy 

i=i 

= (^)( x )+ y {U(x-y)(r|;wi,xi;...;w n ,x n )(y) 

- (U(x-y)r)(y)}dy. (3.10) 

Since x G v\, . . . , v n in the n-th line of the system can take the values of the inclusions v\, . . . , v n , 
the n-th line actually contains n equations. 



4 Background of analytical micromechanics 

4-1 Approximate effective field hypothesis 

In the current section 4, only statistically homogeneous media subjected to homogeneous 
boundary conditions (2.7) are considered. In order to simplify the exact system (3.10) we now 
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apply the so-called effective field hypothesis which is the main approximate hypothesis of many 
micromechanical methods: 

Hypothesis la, Hla. Each inclusion Vi is located in the field (3.92) 

£i(y) = e(xi) (y G Vi ), (4.1) 
which is homogeneous over the inclusion V{. 

In some methods (such as, e.g., the MEFM) this basic hypothesis Hla is complimented by 
a satellite hypothesis [compare with (3.6)]: 

Hypothesis lb, Hlb. The perturbation introduced by the inclusion v\ at the point y ^ v j 
is defined by the relation 

J U(y - x)y i (x)r(x)dx = v t Ti(y - x.)r,. (4.2) 



Hereafter Tj = (r(x)Vj(x))(j) is an average over the volume of the inclusion (but not over 
the ensemble), ((.))» = (((■))(*))' and ( x *' z G x ^ ^> x i>y G v i) 

Ti(x-x.) = (U(x - z)) (l) , Ty( Xj - x,) = (Ti(y - Xi )) (j) . (4.3) 

If x G then Tj(x — xj = — (Uj) _1 Pj =const., where the tensor Pj is associated with the well- 
known Eshelby tensor by S, = Pjl/ ) . For a homogeneous ellipsoidal inclusion Vi the standard 
assumption (4.1) (see, e.g., [1]) yields the assumption (4.2), otherwise the formula (4.2) defines 
an additional assumption. The tensors Tjj(Xj — Xj) has an analytical representation for spherical 
inclusions of different size in an isotropic matrix (see for references [1]). 

According to hypothesis Hla and in view of the linearity of the problem there exist constant 
fourth and second-rank tensors Aj(x), Rj(x), such that 

e(x) = Aj(x)e(xj), r(x) = R(x)e(xj), x G v i: (4.4) 

where Vi C and R(x) = L^Aj(x). According to Eshelby's [13] theorem there are the follow- 
ing relations between the averaged tensors (4.4) R = UjP^ 1 (I — Aj), where g« = (g(x))(j) (g 
stands for Aj,R). For example, for the homogeneous ellipsoidal domain V{ (2.4) we obtain 

Aj = ^1 + PjL^ . In the general case of coated inclusions v^, the tensors Aj(x) can be found 
by the transformation method by Dvorak and Benveniste [14] (see for references and details [1]). 

4-2 Closing hypothesis 
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For termination of the hierarchy of statistical moment equations (3.10) we will use the 
closing EFH called the "quasi-crystalline" approximation by Lax [15] which in our notations has 
a form 

Hypothesis 2, "quasi-crystalline" approximation. It is supposed that the mean value 
of the effective held at a point x 6 Vi does not depend on the stress held inside surrounding 
heterogeneities Vj ^ vf. 

(WiWlvi^^Vj,^) = (Wi), x G v { . (4.5) 



In the framework of the hypothesis HI, substitution of the solution (4.4) into the first 
equation of the system (3.9) at n = 1 and at the EFH H2 leads to the solution (x G vi) 



<e) i = R- 1 YR(e), 
(r) 1 (x)=R(x)R- 1 YR( £ ), 
L * = L (0) +YRc « 



(4.6) 
(4.7) 
(4.8) 



where the matrix Y determines the action of the surrounding inclusions on the considered one 
and has an inverse matrix Y _1 given by 



(Y 1 )=l-V i nJ T iq (Xi -X q )ip(v q ,Xq\;Vi,Xi) - Tj(Xj 



x„ n 



(i) 



dx n 



(4.9) 



General case of the closing hypothesis taking n interacting heterogeneities is considered in Chap- 
ter 8 in [1]. 



4-3 Hypothesis of "ellipsoidal symmetry" of composite structure 



To make further progress, the hypothesis of "ellipsoidal symmetry" 1 for the distribution of 
inclusions attributed to Willis [5] is widely used: 

Hypothesis 3, "ellipsoidal symmetry". The conditional probability density function 
ip(vj,Xj |;uj,Xj) depends on Xj — Xj only through the combination p = |(a^-) _1 (xj — Xj)|: 

tp(vj,Xj \;vi,Xi) = h(p), p=\ (a?-) -1 ( x j ~ x *) I ( 4 - 10 ) 

where the matrix (a^) -1 (which is symmetric in the indexes i and j, a^- = a°J dehnes the 
ellipsoid excluded volume vfj = {x : |(a? ? )~ 1 x| 2 < 1}. 

For spherical inclusions the relation (4.10) is realized for a statistical isotropy of the compos- 
ite structure. It is reasonable to assume that (a- 1 ) -1 identifies a matrix of affine transformation 
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that transfers the ellipsoid being the "excluded volume" ("correlation hole") into a unit 
sphere and, therefore, the representation of the matrix Y can be simplified: 

(Y- 1 ) = I - c«RP°, (4.11) 

where for the sake of simplicity of the subsequent calculation we will usually assume that the 
shape of "correlation hole" vfj does not depend on the inclusion vf vfj = vf and P^- = P° = 

PK°). 

The concept of the EFH (even if this term is not mentioned) in combination with subsequent 
assumptions (e.g., mentioned above) totally dominates (and creates the fundamental limitations) 
in all four groups of analytical micromechanics in physics and mechanics of heterogeneous me- 
dia: model methods, perturbation methods, self-consistent methods (e.g., Mori-Tanaka, MT, 
approach, and the MEFM), and variational ones (see for references [1]). 



5 Background of computational analytical micromechanics 

5.1 A single inclusion subjected to inhomogeneous prescribed effective field 

In the current subsection we will consider a satellite problem whose solution will be used 
for estimation of effective properties of composites in Subsection 5.2. Namely, let the inclusions 
Vi be fixed and loaded by the inhomogeneous effective field £«(x). Then we used the known 
regularized integral equation 

r(x) =Ti(x) + J Ki(x,y)[r(y) - T (x)]dy, x € v h (5.1) 

where Tj(x) = Ej(x)e(x), (x € Vi) is called the effective stress polarization tensor in the inclusion 
Vi, and (no sum on i) 

K t (x,y) = E,(x)U(x-yMy), (5.2) 
E,(x) = MxMl + P^Mx)]- 1 . (5.3) 

Here the tensor P°(x) = - / F i °(y)U(x - y)dy can be estimated, e.g., by the FEA and assumed 
to be known. 

We formally write the solution of Eq. (5.1) as 

T = d*Ti, (5.4) 

where the inverse operator d = (I — /Q)^ 1 will be constructed by the iteration method based 
on the recursion formula 

r [ fc+1 l =Ti + KiT^, (5.5) 

the convergence of which is analyzed in [1]. Here the integral operator 7Q has the kernel formally 
represented as 

7Q(x, y) = Ki(x, y) - <5(x - y) J V-(z)K,(x, z)dz, (5.6) 
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and one used an initial approximation 

r[°](x)=r i (x), (5.7) 

which is exact for a homogeneous ellipsoidal inclusion subjected to remote homogeneous stress 
field e(x) = s = const. 

The solution (5.4) allows us to state that the linear operators C e and C T describing a 
perturbation of the strain fields inside and outside the inclusion Vi (x € R d ) 

|u(x-y)V-(y)r(y)dy = e(x) -s,(x)^(^)(x)^[(r)(x), (5.8) 

A e te)(x) = /U(x - y)d * (E i e)(y)y j (y)dy, (5.9) 

£[(r)(x) = y U(x - y)r(y)y i (y)dy. (5.10) 

are constructed. 

5.2 Estimation of effective elastic moduli 

The new general integral equation (3.5) can be rewritten in terms of the operator represen- 
tation C T (5.8) 

e(x) = (e)(x) + J [£ t (t)(x) - (^(r))(x)]dy. (5.11) 

For statistically homogeneous media subjected to homogeneous boundary conditions (2.7) 
and in the framework of the quasi-crystalline approximation (4.5), conditional averaging of Eq. 
(5.11) can be solved by the iteration method 

(s)!" +1] (x) = ( £ ) + |^((r)W)(x)[^,x ? |;^,x i )-nW(x ? )]dx g , 

(r)[" +1 ](x) = 7^*<e)[" +1 l(x), (5-12) 

where TLq = L^(I + CV). Generalization of Eq. (5.12) to the cases of both the statistically 
inhomogeneous media and other multiparticle closing assumptions (see [1]) is obvious. A con- 
vergence of the sequence (r^^x) (5.12) is analyzed analogously to the sequence (5.5). An 
initial approximation (r)'°'(x) is defined by the classical approach (4.7) and (4.11). It sug- 
gests the Neumann series form for the solution (e)["'(x) — ► (e)j(x) (as n — > oo) of (5.12) and 
(r)!" ] (x)=A*(E J (s)f ] )(x): 

(r) i (x) lim (rW) ,(x) = R*(x)(e), (5.13) 

which yields the final representations for the effective properties 

L* = + (R*V). (5.14) 
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6 Qualitative comparison of the classical and new approaches 

The hypothesis HI is widely used (explicitly or implicitly) for the majority of the methods 
of micromechanics even if the term "effective field hypothesis" is not indicated. For example, 
Buryachenko [1] demonstrated that hypothesis HI is exploited in the effective medium method, 
generalized self-consistent method, differential methods, Mori-Tanaka method, the MEFM, con- 
ditional moments method, variational methods, and others. These are a lot of other methods 
using the hypothesis HI differ one from another by some additional specific assumptions used 
at the analysis of the initial integral equations either Eqs. (3.5), (3.7), or (3.8). 

The differences of Eqs. (3.5), (3.7), and (3.8) are fundamental for subsequently solving 
the truncated hierarchy (3.10) involving a rearrangement of each appropriate equation before 
it is solved. The most successful rearrangement are those which make the right-hand side of 
the coupled equations reflect the detailed corrections to that basic physics. So, Eq. (3.8) was 
obtained by subtracting the difficult state at infinity from equation (3.1), i.e. roughly speaking 
the constant force-dipole density expressed through an alternative technique of the Green's 
function. This dictates the fundamental limitation of a possible generalization of Eq. (3.8) to 
both the FGMs and inhomogeneous boundary conditions. The mentioned deficiency of Eq. (3.8) 
was resolved by Eq. (3.7) which the renormalizing term provides an absolute convergence of the 
integral in Eq. (3.7) at |x — y | — > oo for the general cases of the FGMs. However, the same term 
in Eq. (3.7) is used in a short-range domain |x — y| < 3a in the vicinity of the point x £ w. A 
fundamental deficiency of Eq. (3.7) is a dependence of the renormalizing term U(x — y)(r)(y) 
[obtained in the framework of the asymptotic approximation (3.6)] only on the statistical average 
(r)(y) while the renormalizing term (U(x — y)r)(y) in Eq. (3.5) explicitly depends on on details 
distribution {r\v q ,x q )(y) (y £ v q ). What seems to be only a formal trick [abandoning the use 
of the approximations (3.6)] is in reality a new background of micromechanics [defining a new 
field of micromechanics called computational analytical micromechanics, CAM] which yields the 
discovery of fundamentally new effects even in the theory of statistically homogeneous media 
subjected to homogeneous boundary conditions. So, the final classical representation of the 
effective properties (4.8) depends only on the average strain concentrator factor Aj while the 
effective properties (5.14) implicitly depend on the inhomogeneous tensor Aj(x). Moreover, the 
detected dependence of the effective properties (5.14) on the detailed strain concentrator factors 
Aj(x) rather than on the average values Aj allows us to abandon the hypothesis Hlb whose 
accuracy is questionable for the noncanonical inclusions. In such a case the statistical average 
effective field estimated by Eq. (5.12) is found to be inhomogeneous that discards the hypothesis 
Hla. Thus, the CAM does not involve the hypotheses Hla, Hlb, and H3 as contrasted to 
the classical analytical micromechanics. Only the closing assumption H2 (or its multiparticle 
generalizations, see [1]) is exploited in CAM. 

It should be mentioned, that the domain of the operator £g((l)j^)(x) (5.12) is a whole space 
x G R d , and, because of this, some points of the area x G v i in Eq. (5.12) can be uncovered 
by the heterogeneities v q and, therefore, the effective strain (e)| n+1 '(x) (5.12) will depend on 
the strain perturbations C q ({e) q n ^)(-x) in the vicinity x G v ® of the area v q rather than only 
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on stress distributions in the inhomogeneity v q . In particular, for well-stirred approximation of 
the binary correlation function <p(v q , x 9 |vj, Xj) of ellipsoidal inclusions, v® is expressed by the 
Minkowski addition v ® = v® © while for the spherical inclusions v® = {x| |x — x g | < 3a}. 
Thus, we obtain a fundamental conclusion that effective moduli (5.14) in general depend not 
only on the strain distribution inside the inhomogeneities but also on the strains in the vicinities 
of heterogeneities. Then the size of the excluded volume as well as the RDF will impact on 
the effective field (5.12) even in the framework of hypothesis H3. Indeed, if the radius of the 
excluded volume v® for the spherical inclusions increases from 2a to 3a then the long distance of 
the influence zone v® of the inhomogeneity v q on the effective field (e)j(x) will increase from the 
value |x — x ? | = 3a to |x — x g | = 4a. A larger difference between the backgrounds (3.5) and (3.8) 
was obtained for composites reinforced by either nonellipsoidal or inhomogeneous inclusions 
demonstrating essentially inhomogeneous stress distribution inside isolated heterogeneities even 
in the framework of the hypothesis HI. 

Quantitative estimations of the analyses presented above are under progress for some par- 
ticular cases of fiber composites and will be considered in other publications. 
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